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Resumen 


Considerando la fuerte variabilidad espacial del fenómeno de la 
precipitación en el estado de Nuevo León, México, se utilizaron las 
técnicas geoestadísticas de kriging ordinario y simulaciones estocásticas, 
con el objetivo de estimar la distribución espacial de la precipitación 
media anual y para las estaciones de mayor precipitación (primavera y 
verano). Los datos para el estudio corresponden a 95 estaciones 
climatológicas, con al menos 25 años de registro en el periodo 
comprendido de 1930 a 2014, se incluyen datos de estaciones vecinas 
de la NOAA y estados colindantes a Nuevo León. El uso sistemático de 
estas técnicas permitió establecer la dependencia espacial de los datos 
de precipitación, por medio de modelos variográficos isotrópicos y 
anisotrópicos, los cuales se aplicaron en las estimaciones y 
simulaciones. Para evaluar la calidad de las estimaciones obtenidas se 
aplicó la técnica de validación cruzada. Los resultados obtenidos son 
consistentes con el comportamiento estacional de la precipitación en las 
regiones climáticas del norte y noreste de México, considerando su 
variabilidad espacial, y la disposición y calidad de la información de las 
estaciones climatológicas en el área de estudio. 


Palabras clave: geoestadística, estimación espacial, precipitación, 
Kriging, simulaciones, anisotropía, Nuevo León. 
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Introducción 


Los campos en los que se requiere conocer la distribución espacial de la 
precipitación son diversos, entre éstos se encuentran la hidrología, 
agricultura y climatología; en hidrología, es un insumo importante para 
el modelado hidrológico, en la predicción de eventos extremos, como 
crecientes en el flujo superficial de agua e inundaciones, la presencia de 
lluvia durante el año es muy importante, ya que en cierta medida influye 
en el sector agrícola (Englehart 8 Douglas, 2000), en especial, la 
agricultura de temporal, ya que es dependiente de la precipitación y por 
tanto, altamente vulnerable a la variabilidad interanual e interestacional 
de la misma. Por otro lado, el cambio climático tiene impactos directos 
sobre el régimen de precipitación, lo cual afecta en la gestión de los 
recursos hídricos, de ahí la importancia de analizar y estudiar el cambio 
en el patrón espacial y temporal de la precipitación. 


Los datos utilizados para este estudio provienen de la red de estaciones 
climatológicas de la Comisión Nacional del Agua (CONAGUA), y de la 
NOAA (The National Oceanic and Atmospheric Administration), para el 
caso del estado de Nuevo León se encuentran distribuidas 
irregularmente, y en algunas de estas estaciones, la información 
disponible es insuficiente para caracterizar la alta variabilidad de la 
precipitación y su distribución espacial. Por lo tanto, es necesario 
estimar la precipitación en áreas donde no se dispone de información, 
utilizando datos de las estaciones vecinas. Varios autores como 
Goovaerts (1999); Mirás-Avalos, Paz-González, Vidal-Vázquez y Sande- 
Fouz (2007); Coulibaly y Becker (2007), han demostrado que los 
métodos geoestadísticos proporcionan mejores estimaciones de la 
precipitación, en comparación con otras técnicas. 


Los objetivos principales de este estudio son: 
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. Analizar y modelar la variabilidad espacial de la precipitación media 
mensual y anual. 

. Construir superficies estimadas y de sus errores, utilizando kriging 
ordinario y simulaciones secuenciales gaussianas para la media anual y 
los meses de mayor precipitación, agrupados en dos periodos 
estacionales en Nuevo León, Primavera (febrero 16.8 mm - mayo 55.6 
mm) y Verano (junio 64.8 mm - septiembre 117.7 mm). 


Evaluar los resultados de las técnicas utilizadas. 


Descripción de los datos y área de estudio 


En general, en México, la distribución de precipitación varía en espacio y 
tiempo; es heterogénea a lo largo del año, geográficamente aumenta en 
dirección norte-sur, debido a la influencia de la latitud. Del mismo modo, 
está gobernada en gran medida por la proximidad al Océano Pacífico y 
Golfo de México (Campos-Aranda, 1998), a la orografía del país y a los 
rasgos de la circulación atmosférica (García, 2003). 


El estado de Nuevo León se localiza en el noreste de México, entre los 
meridianos 98% 26' y 101% 14” de longitud oeste, y entre los paralelos 
23% 11' y 27% 49” de latitud norte (Figura 1), cuenta con una superficie 
de 64 000 km? con clima seco. Son tres las provincias fisiográficas que 
lo atraviesan de manera parcial, la Sierra Madre Oriental, las grandes 
llanuras de Norteamérica y la llanura costera del golfo, el rango de 
altitud oscila entre 70 a 3 500 metros. 
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Zacatecas Leyenda 


Altitud (m) 
[E Menor de 400 
[1] 40 - 1.000 
g [11 1.000 - 1,500 
San Luis Potosí [E 1500 - 2.000 
[E] Vayor de 2,000 


Figura 1. Ubicación del área de estudio y distribución de estaciones climatológicas 
utilizadas para la estimación de la precipitación en el estado de Nuevo León, México. 


Los climas secos y semisecos se distribuyen principalmente en la región 
nororiental, la cual forma gran parte de la Gran Llanura de 
Norteamérica; y en la región suroccidental, separada de la primera, por 
las alturas de la Sierra Madre Oriental. En áreas menores de la región de 
la sierra, en la zona del centro y sur de la entidad, y en gran parte de la 
cuenca del Río San Juan se registran los climas semicálidos, templados y 
semifríos (INEGI, 1986), como se observa en la Figura 2. 
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A (ayciwo) Semicálido subhúmedo, con lluvias en verano, de menor humedad|_—] BS1kw Semiseco templado, con lluvias en verano 
EA (ajciws) Semicálido subhúmedo, con lluvias en verano, humedad meda [7] BS1kxX" Semiseco templado, con lluvias escasas todo el año 
MA (Aca) Semicálido subhúmedo, con lluvias en verano, más humedo | BWhw Muy seco semicálido, con lluvias en verano 
(Ne Semicálido subhúmedo, con lluvias escasas todo el año —] Bwhx Muy seco semicálido, con lluvias escasas todo el año 
BS0(h) Seco muy cálido y cálido, con lluvias en verano | C(EXw1) Semifrío subhúmedo, con lluvias en verano, humedad media 
[EJ esonw Seco semicálido, con lluvias en verano EJ] C(ExXw2) Semifrío subhúmedo, con lluvias en verano, más húmedo 
EY] BSOhx Seco semicálido, con lluvias escasas todo el año MA crex Semifrío subhúmedo, con lluvias escasas todo el año 
BSOkw Seco templado, con lluvias en verano MI oo) Templado subhúmedo, con lluvias en verano, menor humedad 
BS0kx' Seco templado, con lluvias escasas todo el año MN cws) Templado subhúmedo, con lluvias en verano, humedad media 
MES Semiseco muy cálido y cálido MA ca) Templado subhúmedo, con lluvias en verano, mayor humedad 
BS1hw Semiseco semicálido, con lluvias en verano | [ey Templado subhúmedo, con Ilwias escasas todo el año 


Figura 2. Tipos de climas de Nuevo León, México. 


La precipitación en el estado está dominada por situaciones estacionales 
y con una gran variabilidad; por ejemplo el área se ve afectada por 
fuertes lluvias, causadas principalmente por tormentas tropicales y 
huracanes, debido a su cercanía al Golfo de México y la influencia de 
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frentes fríos; en general, la precipitación es escasa, aunque cuenta con 
regiones que registran lluvias anuales mayores a 1 100 mm, la media 
anual para Nuevo León en el periodo de estudio es de 541 mm. 


El periodo de 1930 a 2014 fue seleccionado para la integración de los 
datos de precipitación acumulados mensual y anualmente, considerando 
al menos 25 años de registro, para cada una de las estaciones 
climatológicas utilizadas; 56 se encuentran ubicadas en Nuevo León, y 
con la finalidad de contar con una mayor cobertura espacial se 
incorporaron 11 de Coahuila, 17 de Tamaulipas, 4 de San Luis Potosí, 
una de Zacatecas y 6 de Texas, en Estados Unidos de América; tomando 
en cuenta el criterio de mantener una distancia promedio menor de 50 
km al límite del estado de Nuevo León, obteniendo un total de 95 
estaciones, a partir de las cuales se derivaron las variables PMA, 
Primavera y Verano, de acuerdo con la Figura 1. 


Metodología 


La aplicación de la metodología geoestadística a un conjunto de datos, 
involucra tres pasos: análisis exploratorio de datos, análisis estructural, 
y la estimación y/o simulación. En el análisis exploratorio se aplican 
procesos de control de calidad y criterios para la integración de las 
series de datos climáticos, aplicando las recomendaciones publicadas 
por la WMO (World Meteorological Organization) En los documentos, 
Calculation of monthly and annual 30-year standard normals (WMO, 
1989) y Guide to Climatological Practices (WMO, 2011), se analiza la 
estadística básica de los datos mensuales y anuales integrados. 


El análisis estructural consiste en analizar los datos para conocer su 
variabilidad espacial, evaluar la presencia de anisotropía y determinar 
los modelos variográficos para cada variable, para posteriormente 
proceder a validarlos, mediante la técnica de validación cruzada, la 
tercera etapa es la estimación y simulación, en la cual se obtienen las 
superficies interpoladas. En este proyecto, se aplicó además del kriging 
ordinario, la técnica de simulación secuencial gaussiana, la cual consiste 
en que un nuevo valor simulado que se obtiene a partir de una función 
de distribución de probabilidad, utilizando valores muestreados y valores 
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previamente simulados en una vecindad dada, aplicando una técnica de 
kriging (Chilés 8: Delfiner, 1999). 


Resultados 


Análisis Exploratorio 


El control de calidad aplicado a los datos, consistió en la evaluación de 
los siguientes puntos: 


Consistencia espacial. Se actualizaron las coordenadas de 8 
estaciones visitadas: 5 de Coahuila (5016, 5032, 5035, 5048, 5049) y 3 
de Nuevo León (19012, 19054 y 19069), debido a que presentaban 
inconsistencia en su ubicación. 


Formato. Se depuraron registros repetidos y fechas improbables; por 
ejemplo, en los registros de la NOAA, invariablemente todos los meses, 
de todos los años tienen un formato de 31 días. 


Completitud. Para la integración de los datos anuales, se seleccionaron 
las 95 estaciones que cumplieron con un periodo de al menos 25 años 
de registro, considerando el criterio de excluir los meses con más de 5 
registros diarios faltantes o más de 3 registros diarios continuos. 


La estadística básica para las variables PMA, Primavera y Verano, se 
describe en la Tabla 1. Las tres variables de precipitación tienen una 
distribución asimétrica positiva, por lo que es necesario evaluar si en 
estos casos se presenta el conocido “efecto proporcional”, que es una 
forma particular de heterocedasticidad (la variabilidad de los datos 
cambia a lo largo del área de estudio), en particular para distribuciones 
asimétricas positivas, la varianza local se incrementa conforme se 
aumenta su media local (Goovaerts, 1997), este efecto proporcional 
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puede hacer el variograma experimental no interpretable (Grimes 48 
Pardo-Igúzquiza, 2010). 


Tabla 1. Estadística básica de las variables analizadas. DE (desviación estándar), CV 
(coeficiente de variación). 


Media DE cv Min Qi Mediana Q3 Max 
Variable Asimetría | Curtosis 
(mm) (mm) (mm) | (mm) (mm) (mm) (mm) (mm) 

PMA 539.2 | 199.8 | 0.4 | 212.7|394.5| 507.6 | 661.9 | 1,110 0.77 3.2 
log PMA 6.2 0.4 0.1 5.4 6.0 6.2 6.5 7.0 -0.1 2.7 
Verano 308.3 | 126.6 | 0.4 | 114.8 |222.0| 285.6 | 367.1 | 700.9 0.9 3.7 

log 
O 0.4 | 0.1 | 4.7 | 5.4 57 5.9 6.6 -0.05 2.6 

Primavera | 123.5 | 44.8 0.4 40.6 98.5 117.7 | 147.7 | 243.1 0.5 2.9 


Para conocer cómo están relacionados los cambios en la varianza local, 
con los cambios en la media local, es común utilizar gráficas de 
dispersión calculadas a partir de ventanas móviles (Goovaerts, 1997). 
Se calcularon las estadísticas para ventanas móviles con un área de 100 
km?, con sobreposición de 25 km”; se seleccionaron 18, que incluyen al 
menos 10 datos. La sobreposición de las ventanas implica que algunos 
datos fueron utilizados en repetidas ocasiones para el cálculo de las 
medias y varianzas locales, esto no tiene relevancia desde el momento 
en que sólo se quiere establecer si existe o no el efecto proporcional 
(Isaaks €: Srivastava, 1989), el resultado obtenido para PMA y Verano 
muestra que, las varianzas y medias locales tienen baja correlación 
(Pranx = 0.31) por lo que se puede asumir la homocedasticidad, en tanto 
que Primavera presenta un efecto proporcional, debido a una correlación 
mayor (Pranx = 0.62); en la Figura 3 se muestran las gráficas de 
dispersión. 
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8] Prank = 0.31 
400 500 e 700 BR PMA 
Prank = 0.31 E Ñ Prank = 0.62 
Verano Primavera 
Figura 3. Gráficas de dispersión para 18 ventanas móviles con al menos 10 datos por 


ventana. 


Se sabe que las estaciones climatológicas son ubicadas en lugares 
donde se cuente con ciertas características para su operación y registro, 
por lo que existen áreas con mayor concentración (clustering), en 
particular cuando se combinan el efecto proporcional y el agrupamiento 
de los datos de la muestra, conducen a serios problemas en la 
interpretación de los variogramas experimentales, una manera de 
conocer el agrupamiento y el efecto proporcional es a través de graficar 
las medias y varianzas locales en función de la distancia (Goovaerts, 
1997). 


Para el caso de las tres variables se observa un comportamiento que 
fluctúa alrededor de un valor unitario, concluyendo que las varianzas y 
las medias no se afectan de manera significativa con la distancia, por lo 
que no se requiere realizar algún proceso para considerar el efecto de 
agrupamiento; por ejemplo, usar solo datos regularmente espaciados 
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(Goovaerts, 1997). En la gráfica de la Figura 4, sólo se muestran los 
resultados obtenidos para PMA. 
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Figura 4. Medias y varianzas como una función de la distancia (lags). Ambos 
estadísticos son normalizados por sus correspondientes media y varianza global. 


Análisis estructural 


Considerando que el grado de asimetría y los datos atípicos afectan de 
manera directa en el modelado de los variogramas, se aplicó una 
transformación logarítmica para disminuir el grado de asimetría y 
analizar el comportamiento de los datos atípicos. La transformación 
logarítmica utilizada para Primavera no mejoró la distribución, por lo 
cual no se consideró su aplicación. Las gráficas de la Figura 5, muestran 
los variogramas experimentales para las variables con asimetría 
moderada (PMA y Verano), se incluyó su transformación logarítmica en 
la misma gráfica, después de reescalar la varianza. Los variogramas 
logarítmicos parecen menos erráticos, con menor nugget para PMA, por 
lo tanto, esta variable se modela sin transformación, mientras que 
Verano sí justifica la transformación, debido a la notable diferencia entre 
los variogramas. 
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Figura 5. Variogramas experimentales omnidireccionales para PMA y Verano, el 
variograma de los logaritmos es reescalado a la varianza de los datos originales. 


El análisis de anisotropía se realizó inspeccionando los variogramas 
experimentales en diferentes direcciones: 0% (N-S), 45(NE-SO), 90%(E- 
O) y 135%(NO-SE) con una tolerancia angular de + 22.5 2 
(Wackernagel, 2003). Como en el área de estudio se encuentra ubicada 
una porción de la Sierra Madre Oriental, la cual influye en los patrones 
de la precipitación, se construyeron variogramas en la dirección NO 20%, 
rumbo asociado a la estructura regional de la sierra citada, así como la 
dirección perpendicular a ésta (NE 70%), y adicionalmente se construyó 
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el variograma omnidireccional (dirección 0%, con una tolerancia angular 
de 909), con la ayuda del software gstat (Pebesma € Wesseling, 1998) y 
el RGEOESTAD (Díaz, Hernández, 8 Méndez, 2012). 


De acuerdo con las gráficas de la Figura 6, a partir de los puntos de 
inflexión de los variogramas experimentales, se aprecian dos escalas de 
variación espacial en la dirección NE 70% (perpendicular a la estructura 
de la Sierra Madre Oriental): una con un rango de 75 km y otra de 125 
km aproximadamente, para una distancia mayor a 125 km, la estructura 
de los variogramas se vuelven erráticas en esa dirección, debido a la 
disminución del número de pares que contribuyen a los valores de las 
varianzas, para el caso de la dirección NE 20%, los variogramas 
muestran una mayor continuidad espacial (menor varianza) y siguen un 
patrón similar al del variograma omnidireccional. 
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Figura 6. Variogramas PMA, Primavera y Log Verano para las direcciones NE 70%, NO 

20% con tolerancia angular de 22.5 líneas discontinuas, el variograma omnidireccional 

(09) con tolerancia angular de 90% en línea continua. Las etiquetas numéricas detallan 
el número de pares. 
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Otra manera de evaluar la anisotropía es a través de la elaboración de lo 
que se conoce como mapa variográfico (Isaaks 8. Srivastava, 1989), en 
el cual se grafican los valores del variograma experimental y el centro 
del mapa corresponde al origen del variograma, en la Figura 7 se 
observa una dirección preferencial hacia el NO-SE para Log Verano y 
PMA, aproximadamente alineada a la que presenta la Sierra Madre 
Oriental; para el caso de Primavera se aprecia en menor medida esta 
dirección preferencial. 


T T T T Y T T 
300000  -200000  -100000 0 100000 200000 300000 


Log Verano 
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Primavera 
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Figura 7. Mapas variográficos, PMA, Log Verano y Primavera, valores 
estandarizados por la varianza de la muestra. 


Tomando en cuenta los resultados arrojados por el mapa variográfico, 
se utilizaron los modelos anisotrópicos para las variables PMA y log 
Verano. 


Una vez definidas las direcciones de los variogramas se ajustaron 
modelos teóricos esféricos, con la técnica de mínimos cuadrados, con el 
apoyo del software RGEOESTAD (Díaz et a/., 2012) y gstat (Pebesma € 
Wesseling, 1998), se obtuvieron los parámetros variográficos: nugget, 
sill y rango; en la Tabla 2 se muestran los datos obtenidos para las tres 
variables analizadas, y la Figura 8 ilustra los modelos esféricos 
ajustados, incluyendo el omnidireccional. Los modelos ajustados para la 
dirección NE 70% no se muestran en este documento. 


Tabla 2. Parámetros variográficos NO 20%, NE 70% y omnidireccional, modelados para 
las variables Primavera, PMA y Verano logarítmico. 


Dirección NO 20? NE 702 Omnidireccional 
Variable 
PMA Primavera Log PMA Primavera Log PMA Primavera Log 
Verano Verano Verano 
A | 6500 310 0.04 | 16860 | 537 0 9 570 398 0.03 
Sill 
41 118 2140 0.18 47 360 2311 0.17 33 850 1752 0.15 
(mm?) 
Rango 
150 215 230 62 90 44 77.2 103.6 75.7 
(km) 
Nugget/ 
Sil 0.16 0.14 0.22 0.36 0.23 0 0.28 0.23 0.2 
] 
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Figura 8. Variogramas direccionales ajustados para PMA, Primavera y Log 
Verano. a) Dirección NO 20%; b) Omnidireccional. 
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Todos los variogramas para la dirección NO 20%, exhiben un 
comportamiento lineal en el origen, así como pequeños efectos nugget, 
cuya proporción con respecto al sill se encuentra en el intervalo del 16% 
al 23%, y representan una discontinuidad en el origen que es atribuible 
a errores de medición y a la variación a distancias inferiores del 
intervalo de muestreo (28.4 km). 


Los variogramas en la dirección NE 70% se modelaron para obtener el 
factor de anisotropía A, definido como la razón del rango menor (NE 
70%) y el rango mayor (NO 209). 


Validación cruzada 


La técnica de validación cruzada consiste en sacar una observación y 
estimar su valor con el resto de las observaciones, esto se realiza 
repetidamente hasta concluir con todas las observaciones. 


Los errores son la diferencia entre los valores observados y estimados, y 
se calculan por medio de la siguiente expresión: 


E=2Z06)£ i=Linñn (1) 


Donde e son los errores; Z el valor observado y Z* el estimado. 


Existen una serie de indicadores que miden la calidad de la estimación 
de manera global, en este trabajo se consideran los siguientes: 


1 
ME a i=1 € (2) 


RMSE = [Ly e? (3) 
n 


La Tabla 3 muestra los resultados de la validación. 


Tabla 3. Resultados de los errores obtenidos de la validación cruzada. 
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Variable Modelo 


anisotrópico 


isotrópico 


anisotrópico 


Verano 


isotrópico 


anisotrópico 
Primavera 
isotrópico 


De acuerdo con el RMSE obtenido, el modelo anisotrópico mejoró la 
estimación de PMA y Primavera, para Verano resultó mejor la estimación 
del modelo isotrópico, de acuerdo con los resultados obtenidos en la 
validación cruzada. 


Métodos de estimación y simulación 


Los mapas de precipitación estimada y error de la estimación mostrados 
en la Figura 9, Figura 10 y Figura 11 se realizaron con el software 
ArcMap, utilizando la extensión geostatistical analyst (ESRI, 2016), 
utilizando los modelos variográficos ajustados, aplicando las técnicas de 
kriging ordinario y simulaciones secuenciales gaussianas, estas técnicas 
son discutidas y analizadas ampliamente en los trabajos de algunos 
autores como Chiles y Delfiner (1999), Goovaerts (1997), Wackernagel 
(2003), e Isaaks y Srivastava (1989). 
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Coahuila de Coahuila de Coahuila de 
Zaragoza Zaragoza Zaragoza 
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San Luis Potosí San Luis Potosí Km San Luis Potosí 


MD 2104-2025 
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ENNCIEI 


(d) Error modelo isotrópico (e) Error modelo anisotrópico 


San Luis Potosí 


Figura 9. Mapa de Nuevo León de las estimaciones, simulación y error de la 
estimación para PMA, en el periodo 1930-2014. 
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Zaragoza Zaragoza 
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(c) simulación 


Coahuila de 
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San Luis Potosí San Luis Potosi 


(d) Error modelo isotrópico (e) Error modelo anisotrópico 


Figura 10. Mapa de las estimaciones de Nuevo León, simulación y error de la 
estimación para Primavera, en el periodo 1930-2014. 
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(b) Modelo anisotrópico 
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(d) Error modelo isotrópico (e) Error modelo anisotrópico 
Figura 11. Mapa de las estimaciones de Nuevo León, simulación y error de la 


estimación para Verano, en el periodo 1930-2014. 
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Para la aplicación de la técnica de simulación secuencial gaussiana, se 
utilizaron los modelos variográficos isotrópicos, aplicando a los datos 
una transformación de anamorfosis. Se utilizó el programa de código 
abierto gstat (Pebesma € Wesseling, 1998), para realizar 100 
simulaciones para cada una de las variables y con el promedio de estas 
simulaciones se generaron los mapas, utilizando la extensión 
geostatistical analyst de ArcMap (ESRI, 2016). 


En el caso de PMA, el modelo anisotrópico, reprodujo de mejor 
manera la distribución espacial de los datos de precipitación, 
enfatizándose la dirección preferencial a la Sierra Madre Oriental, en la 
parte norte del estado se tienen patrones muy similares entre los tres 
modelos, y como es de esperarse, para el caso de las simulaciones, 
produjo fielmente los valores extremos de los datos. 


Conclusiones 


El uso de técnicas geoestadísticas, permitió establecer modelos 
variográficos isotrópicos y anisotrópicos de dependencia espacial para 
los datos de precipitación, los cuales fueron fundamentales para la 
aplicación del kriging ordinario y de simulaciones secuenciales 
gaussianas. 


Las simulaciones secuenciales gaussianas definieron, de manera más 
detallada, los patrones de precipitación estacional para el estado, debido 
a eventos ciclónicos provenientes del Mar Caribe y Golfo de México, con 
una dirección dominante NO-SE y con valores superiores a los 40 mm, 
estos valores contrastan con los resultados de precipitación media 
mensual, interpolados con técnicas de regresión para la región 
hidrológica 24 Río Bravo-Conchos, la cual abarca una porción del estado 
de Nuevo León (Núñez-López et al., 2013), encontrando valores altos 
para los meses lluviosos, predominantemente en las partes altas de la 
Sierra Madre Oriental; sin embargo, no se logra identificar con claridad 
los valores altos de precipitación en las proximidades del área 
metropolitana de Monterrey, pertenecientes a la porción territorial de la 
provincia fisiográfica Llanura Costera Golfo Norte, tal vez debido a la 
cantidad y/o distribución de las estaciones climatológicas utilizadas. 
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Por otro lado, con relación a los valores máximos de precipitación media 
anual obtenidos destaca la región citrícola comprendida por los 
municipios de Allende, Montemorelos, Hualahuises, General Terán y 
Linares, esta porción territorial corresponde a la provincia fisiográfica 
Llanura Costera Golfo Norte, además del Municipio de Santiago y 
perteneciente a la provincia de la Sierra Madre Oriental, cuyos valores 
oscilan entre 800 y 1 110 mm, acordes con el rango de 800 a 1 200 
mm, reportado en Vidal (2005). 


Considerando los valores obtenidos de RMSE, en la validación cruzada, 
el uso de los modelos anisotrópicos mejoraron las estimaciones para 
PMA y Primavera, mientras que para Verano, el modelo isotrópico 
resultó mejor, en todos los casos es poco significativa la mejoría, por lo 
que se consideró usar el modelo isotrópico para las simulaciones. 

A pesar de que el tiempo de procesamiento es mayor en las 
simulaciones secuenciales gaussianas se recomienda su uso, debido a 
que se definieron mejor los patrones de precipitación, que con la 
utilización de kriging ordinario. 

La incorporación de estaciones climatológicas de estados colindantes, 
incluyendo las de la NOAA en Texas, contribuyeron a una mayor y mejor 
cobertura espacial, lo que permitió contar con información suficiente 
para realizar las estimaciones en todo el territorio de Nuevo León. 
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